Search for coherent charged pion production in neutrino-carbon interactions 
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We report the result from a search for charged-current coherent pion production induced by muon 
neutrinos with a mean energy of 1.3 GeV. The data are collected with a fully active scintillator 
detector in the K2K long-baseline neutrino oscillation experiment. No evidence for coherent pion 
production is observed and an upper limit of 0.60 x 10 -2 is set on the cross section ratio of coherent 
pion production to the total charged-current interaction at 90% confidence level. This is the first 
experimental limit for coherent charged pion production in the energy region of a few GeV. 

PACS numbers: 13.15.+g,25.30.Pt,95.55.Vj 



The charged-current (CC) coherent pion production in 
neutrino-nucleus scattering, + A — ► \i~ + 7r + + A, is 
a process in which the neutrino scatters coherently off 
the entire nucleus with a small energy transfer. Such 
a process has been measured in a number of experi- 
ments P, 0, El 0] , providing a test of the partially con- 
served axial- vector current (PCAC) hypothesis [5J. The 
existing data agree with the Rein and Sehgal model 
based on the PCAC hypothesis for neutrino energies from 
7 to 100 GeV, while there exists no measurement at lower 
energies. 

The recent discovery of neutrino oscillations has re- 
newed interest in neutrino-nucleus interactions in the 
sub- to few GeV region. The KEK to Kamioka (K2K) 
long-baseline neutrino oscillation experiment has re- 
ported Q a significant deficit in the forward scatter- 
ing events, which limits the prediction accuracy of the 
neutrino energy spectrum at the far detector. CC co- 
herent pion production is one of the candidate interac- 
tions responsible for this deficit and its study is nec- 
essary to improve the accuracy of the current and fu- 
ture atmospheric/accelerator-based neutrino oscillation 
experiments, which are expected to achieve much im- 
proved statistical precision using interactions of neutrinos 
in the same energy region as K2K. 

This letter presents the result from a search for CC 
coherent pion production by neutrinos in the K2K ex- 
periment. We compare our result specifically with the 
Rein and Sehgal model [|| because it is the only model 
that provides the kinematics of pions and is commonly 
used in neutrino oscillation experiments. 

In the K2K experiment, protons are extracted from the 
KEK 12 GeV proton synchrotron and hit an aluminum 
target. Positively charged secondary particles, mainly pi- 
ons, are focused by a magnetic horn system and decay 
to produce an almost pure (98%) beam with a mean 
energy of 1.3 GeV [8j. The neutrino beam energy spec- 
trum and spatial profile are measured using a set of near 
neutrino detectors located 300 m downstream from the 
proton target. The estimated absolute flux has a large 
uncertainty due to difficulties in the absolute estimation 
of the primary proton beam intensity, the proton tar- 
geting efficiency, and hadron production cross sections. 
Therefore, the ratio of the CC coherent pion to the total 
CC cross section is measured, rather than the absolute 



CC coherent pion cross section. The data used for this 
analysis were collected with one of the near detectors, 
the fully active scintillator detector (SciBar), from Octo- 
ber 2003 to February 2004, corresponding to 1.7 x 10 19 
protons on target (POT). 

The SciBar detector 9] consists of 14,848 extruded 
plastic scintillator strips read out by wavelength-shifting 
fibers and multi- anode photomultipliers. The scintilla- 
tor also acts as the neutrino interaction target; it is 
a fully active detector and has high efficiency for low 
momentum particles. Scintillator strips with dimen- 
sions of 1.3 x 2.5 x 300 cm 3 are arranged in 64 lay- 
ers. Each layer consists of two planes to measure hor- 
izontal and vertical position. The total size of the de- 
tector is 3.0 x 3.0 x 1.7 m 3 , while an inner volume of 
2.6 x 2.6 x 1.35 m 3 (9.38 tons) is used as the fiducial 
volume to reject incoming particles and obtain a flat effi- 
ciency for CC interactions. The minimum reconstructible 
track length is 8 cm. A track finding efficiency of more 
than 99% is achieved for single tracks with a length of 
more than 10 cm. The track finding efficiency for a sec- 
ond, shorter track is lower than that for single tracks due 
to overlap with the first track. This efficiency increases 
with the length of the second track and reaches 90% at 
a track length of 30 cm. 

The NEUT Monte Carlo (MC) simulation program li- 
brary |lfj| is used to simulate neutrino-nucleus interac- 
tions. The CC coherent pion production is incorporated 
in the simulation based on the Rein and Sehgal model , 
which predicts the cross section averaged over the K2K 
neutrino energy spectrum of 2.85 x 10 _40 cm 2 /nucleon for 
carbon. The Llewellyn Smith model 0] and the Rein 
and Sehgal model ^3] are employed for quasi-elastic (QE) 
scattering (v^ + n — > /i~ + p) and CC single pion (l7r) 
production [y^ + N — > /i + N + n), where N is a nu- 
cleon, respectively. The axial vector mass of the nu- 
cleon form factor is set to be 1.1 GeV/c 2 for both QE 
and CC17T interactions For deep inelastic scatter- 

ing (DIS), we use GRV94 nucleon structure functions Q 
with a correction by Bodek and Yang 0] . Nuclear effects 
are taken into account; for the pions originating from 
neutrino interactions, absorption, elastic scattering, and 
charge exchange inside the target nucleus are simulated. 
Pion cross sections are calculated using the model by Sal- 
cedo et al. which agrees well with past experimental 
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data 01- Pi° n interactions outside the target nucleus 
are simulated based on other experimental data |l8|. 

For the present analysis, the experimental signatures 
of CC coherent pion production are the existence of ex- 
actly two tracks, both consistent with minimum ioniz- 
ing particles, and small momentum transfer defined as 
q 2 = (P M — P„) 2 , where and P„ are the four momenta 
of the muon and the neutrino, respectively. According 
to the MC simulation, the dominant background is the 
CC17T production, where the proton is below threshold 
or the neutron is invisible. 

Charged current (CC) candidate events are selected by 
requiring that at least one reconstructed track starting 
in the fiducial volume is matched with a track or hits in 
the muon range detector (MRD) [lj| located just behind 
SciBar (SciBar-MRD sample). This criterion imposes a 
threshold for muon momentum of 450 MeV/c. Ac- 
cording to the MC simulation, 98% of the events selected 
by this requirement are CC induced events, and the rest 
are neutral current (NC) interactions accompanied by a 
charged pion or proton which penetrates into the MRD. 
The contribution from v e is negligible (< 0.4%). The 
momentum of the muon is reconstructed from its range 
through SciBar and MRD. The resolutions for p M and the 
angle with respect to the neutrino beam direction (9^) are 
determined to be 80 MeV/c and 1.6 degrees, respectively. 

From the SciBar-MRD sample, events with two recon- 
structed tracks are selected. The QE candidate events 
are rejected by using kinematic information Q. Events 
in which the shorter track is identified as proton-like 
based on dE/dx information (non-QE-proton sample) are 
also rejected to select the non-QE-pion sample, which 
includes the signal candidates. The particle identifica- 
tion capability is verified using cosmic ray muons and 
the shorter tracks in the QE sample, where the latter 
provides a proton sample with more than 90%. The prob- 
ability to mis- identify a muon track as proton- like is 1.7% 
with a corresponding proton selection efficiency of 90%. 

The CC coherent pion candidates are extracted from 
the non-QE-pion sample. The background events are 
suppressed by requiring that the pion-like track goes for- 
ward. Even if the additional particles in the background 
process are not reconstructed as tracks, they can be de- 
tected as a large energy deposit or additional hits around 
the vertex. Figure ^a) shows a distribution of energy 
deposited in the vertex strip (-Evtx) for the non-QE-pion 
sample. The MC prediction for -Evtx is verified with the 
QE sample, which has no contribution from non-visible 
particles, as shown in Fig. 1(b). We require the events to 
have E vtx less than 7 MeV and no additional hits around 
the vertex strip. 

The value of q 2 reconstructed from p M and M under 
the assumption of QE interaction is denoted q 2 cc , and is 
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FIG. 1: Distribution of -Evtx for (a) non-QE-pion sample and 
(b) QE sample. Black circles: observed data, Histograms: 
MC expectation with breakdown of interaction modes. The 
statistical x 2 /DOF in the selected region of (a), indicated by 
a vertical line, is 30.1/7 (9.8/7) with (without) CC coherent 
pion production. 



calculated using 



Vv 



1 (M 2 - ml) + 2E^(M n -V)- (M n - Vf 

2 -En +(M n -V)+ Pt , cos 9^ 



where M p r n -\ is the proton (neutron) mass, is the 
muon mass and V is the nuclear potential set to 27 MeV. 
The q 2 oc for coherent pion production events, which is ex- 
pected to be very small due to the small scattering angle 
for muons, is shifted from the true q 2 by 0.008 (GeV/c) 2 
with a resolution of 0.014 (GeV/c) . Events are required 
to have a reconstructed q 2 of less than 0.10 (GeV/ c) . 

The background contamination in the final sample is 
estimated by the MC simulation. In order to constrain 
the uncertainties, the q 2 ec distributions of the data in the 
region q 2 cc > 0.10(GeV/c) 2 are fitted with MC expecta- 
tions. The one track sample is used as well as two-track 
QE, non-QE-proton and non-QE-pion samples, and these 
four samples are fitted simultaneously. In the fit, the 
non-QE to QE relative cross section ratio, the magni- 
tude of the nuclear effects and the momentum scale for 
muons are treated as free parameters. Figure [5] shows 
the q 2 cc distributions of the data with the MC simula- 
tion after the fitting. The % 2 value in the regions with 
g r 2 ec > 0.10(GeV/c) 2 at the best fit is 73.2 for 82 degrees 
of freedom (DOF). 

Figure shows the <y 2 cc distribution for the final CC 
coherent pion sample. The number of events in each se- 
lection step is summarized in Table [i] together with the 
signal efficiency and purity. In the signal region, 113 
coherent pion candidates are found. The neutrino en- 
ergy spectra for coherent pion events and the efficiency 
as a function of neutrino energy, estimated using the MC 
simulation, are shown in Fig.0{a) and^Jc), respectively. 
The total efficiency is 21.1%. The expected number of 
background events in the signal region is 111.4. After 
subtracting the background and correcting for the effi- 
ciency, the number of coherent pion events is measured 
to be 7.64 ± 50.40 (stat.), while 470 events are expected 
from the MC simulation. Hence, no evidence of coherent 
pion production is found in the present data set. 
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FIG. 2: The g 2 ec distributions for the (a) ltrack, (b) QE, (c) 
non-QE-proton, and (d) non-QE-pion samples. The statisti- 
cal x 2 / D OF in the region g 2 oc < 0.10 (GeV/c) 2 of (c) and 
(d) are 7.2/2 (2.7/2) and 32.3/2 (1.2/2) with (without) CC 
coherent pion production. 
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FIG. 3: The reconstructed q distribution in the final sample. 







(%) 


(%) 


SciBar-MRD 


10049 


77.9 


3.6 


Two track 


3396 


35.5 


5.1 


Non-QE pion 


843 


27.7 


14.8 


Second track direction 


773 


27.3 


15.8 


No activity around the vertex 


297 


23.9 


28.2 


q? cc < 0.10(GeV/c) 2 


113 


21.1 


47.1 



TABLE I: The number of events, the MC efficiency and purity 
of coherent pion events after each selection step. 
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FIG. 4: Top: The neutrino energy spectra for (a) the coherent 
pion and (b) total CC events. The hatched histograms show 
the selected events. Bottom: The efficiencies as a function of 
neutrino energy for (c) the coherent pion and (d) total CC 
events. All of them are estimated by the MC simulation. 



The total number of CC interactions is estimated by 
using the SciBar-MRD sample. As shown in Table I, 
10049 events fall into this category. Based on the MC 
simulation, the selection efficiency and purity for CC in- 
teractions in the sample are estimated to be 56.9% and 
98.0%, respectively. The expected neutrino energy spec- 
tra and the energy dependence of the selection efficiency 
for CC events are shown in Fig. Efb) and^fd), respec- 
tively. The total number of CC events is obtained to be 
(1.73 ± 0.02(stat.)) x 10 4 . We derive the cross section 
ratio of CC coherent pion production to the total CC 
interaction to be (0.04 ± 0.29(stat.)) x 10~ 2 . 

Systematic uncertainties for the cross section ratio are 
summarized in Tabic ITT1 The major contributions come 
from uncertainties of nuclear effects and the neutrino in- 
teraction models. The uncertainty due to nuclear effects 
is estimated by varying the cross sections of pion ab- 
sorption and elastic scatterin g by ±30% based on the 
accuracy of the reference data [l7|- The uncertainties in 
QE and CC17T interactions are estimated by changing the 
axial vector mass by ± 0.10 GeV/c 2 Q. For DIS, the 
effect of the Bodek and Yang correction is evaluated by 
changing the amount of correction by ±30%. The q 2 ec 



distribution of the non-QE-proton sample (Fig. [^c)) in- 
dicates an additional deficit of background events in the 
region q^ cc < 0.10 (GeV/c) 2 . CC17T interaction domi- 
nates events in this region; its cross section has signif- 
icant uncertainty due to nuclear effects. We estimate 
the amount of possible deficit in the same manner as de- 
scribed in |j| with the one track, QE and non-QE-proton 
samples. We find that a 20% suppression of CC17T events 
for g 2 ruc < 0.10 (GeV/c) 2 is allowed, which varies the 
cross section ratio by±0.14xl0~ 2 . This variation is con- 
servatively treated as a systematic uncertainty. We also 
consider the uncertainties of the event selection, where 
the dominant error comes from track counting, detector 
response such as scintillator quenching, and neutrino en- 
ergy spectrum shape. The total systematic uncertainty 
on the cross section ratio amounts to +0.32/— 0.35 x 10 -2 . 

Our result is consistent with the non-existence of CC 
coherent pion production at K2K neutrino beam energies, 
and hence we set an upper limit on the cross section ratio 
at 90% C.L. : 

cr(CC coherent ir)/cr(^CC) < 0.60 x 10~ 2 . 
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Error source 


Uncertainty of a 


ratio (xl0~ 2 ) 


Nuclear effects 


+0.23 


-0.24 


Interaction model 


+0.10 


-0.09 


CCi7r suppression 


+0.14 




Event selection 


+0.11 


-0.17 


Detector response 


+0.09 


-0.16 


Energy spectrum 


+0.03 


-0.03 


Total 


+0.32 


-0.35 



TABLE II: The summary of systematic uncertainties in the 
(CC coherent pion) / (total CC interaction) cross section ratio. 



For reference, the total CC cross section is calculated as 
1.07 x 10 _38 cm 2 /nucleon in the neutrino MC simulation 
by averaging over K2K neutrino beam energies. 

The obtained upper limit is inconsistent with the 
model prediction by Rein and Sehgal at the level of 2.5 
standard deviations. We assign a 35 % uncertainty to the 
theoretical prediction as described in jfj] . In addition, a fi- 
nite cross section was reported by Aachen-Padova group 
for NC coherent pion production with 2 GeV average 
neutrino energy and with aluminum target |2fJ |. If we 
assume an A 1 / 3 dependence of the cross section (a) and 
er(CC) = 2er(NC) according to the model of Rein and 
Sehgal, the discrepancy between the extrapolation from 
the NC measurement and the present result is as large 
as 3 standard deviations. There are other models pre- 
dicting lower cross sections 0, 153, Hj| , but they do not 
provide the kinematics of pions and it is difficult to test 
them directly. Further theoretical work is necessary to 
construct interaction models which explain these exper- 
imental results. The non-existence of CC coherent pion 
production has given a solution to the low-g 2 discrepancy 
observed in K2K. It also reduces the uncertainty on the 
cross section in the relevant q 2 region, which is crucial 
for the future neutrino oscillation experiments. 

In summary, we report on a search for CC coherent 
pion production by muon neutrinos with a mean energy 
of 1.3 GeV. The data analyzed correspond to 1.7 x 10 19 
POT recorded with the K2K-SciBar detector. No ev- 
idence of CC coherent pion production is found and an 
upper limit on the cross section ratio of CC coherent pion 
production to the total CC interaction is derived to be 
0.60 x 10 -2 at 90% C.L. This result is the first experimen- 
tal limit for CC coherent pion production by neutrinos 
with energies of a few GeV. 
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